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Abstract. CMB anisotropics are modified by the weak lensing effect of intervening 
large scale structures on the photon path from the last scattering surface to the 
observer. This has to be accounted for when observational data of sensitive experiments 
are used to constrain cosmological models. A common approximation to analyze 
the CMB angular power spectra is to include only the Gaussian part of the lensing 
correction and to ignore the non-gaussian terms in the error covariance matrix of 
the spectra. In order to investigate the validity of this approximation, we computed 
these non-Gaussian terms by using a perturbative expansion method. We present a 
graphical method to write down any N-point correlation functions at any order in 
lensing. We use a pedagogical approach to demonstrate that neglecting non-gaussian 
terms is an accurate approximation for all polarizations but B, and it will remain so 
even for the analysis of very sensitive post-Planck experiments. For the B polarization, 
non-gaussian contributions up to order 4 must be taken into account. 



PACS numbers: 98.80.-k,98.70.Vc,98.62.Sb,98.65.Dx 



Probing inflation with CM B polarization : weak lensing effect on the covariance of CMB spectra2 
1. Introduction 

Measurements of the temperature anisotropies and polarization of the Cosmic Microwave 
Background (CMB) provide very valuable limits and constraints on our models of 
the early universe pQ. This will be even more so with the increased precisions of 
measurements that future CMB experiments promis^il [2]. However, while the hopes of 
high precision cosmology become true, second order effects that were formerly neglected 
or ignored have now to be taken into account. Among them, one of the most important, 
in particular at small angular scale, is the gravitational shear effect of the large scale 
structure, also commonly called the weak lensing effect [3J. This effect generates B modes 
[I], which are of great interest for testing theories of the early universe. Indeed, scalar 
primordial perturbations do not generate B modes, whose detection could therefore 
be considered as the "smoking gun" of the primordial gravitational wave background 
generically predicted by inflation theories j5]. In addition to determining the energy 
scale of inflation, a detection of the primordial B-mode would also test other aspects 
of the early universe physics, for example the presence of cosmic strings [6j [7]. From 
that point of view, the lensing effect is an annoying foreground limiting the detection 
capability of primordial B modes. It is however quite interesting in its own right, since 
it allows reconstructing the (projected) matter power spectrum as well as extracting 
information on various aspects of high energy physics such that neutrino masses [8j E] . 

The gravitational deviations of photons along their paths as they cross the potential 
wells of the large scale structures has been studied in great details in the context of 
CMB. The gravitational shear shifts power between scales and creates mode couplings 
in the CMB power. It also deviates the distribution of temperature and polarization 
anisotropies from the Gaussian statistic [TO] . The specific and predictable signature of 
the lensing effect has been used to propose various ways of detecting and reconstructing 
its contribution to the temperature and polarization of the CMB (TT], [12], [T3J [HI [15] . 

Until recently however, few studies had been devoted to a detailed assessment 
of the impact of weak lensing on the cosmological parameter estimations from CMB 
data. Indeed, analyzes had in general been concerned with including lensing corrections 
to the power spectra of CMB temperature and polarization anisotropies, but they 
neglected the non-gaussian corrections. Recently A. Lewis showed numerically that this 
approximation is valid at the sensitivity level of the Planck mission [16]. This paper 
presents an analysis which demonstrates that the lensing induced non-gaussian terms of 
the error bars on the CMB anisotropies power spectrum can be neglected as compared 
to the Gaussian sample variance terms except for the B-mode power spectrum. While 
developed independently, this works follows the line of the recently published [T7|, [18] 
and fully confirm their relevant results. Given the availability of the numerical results 
of these references, we rather focus here on the physical origin of the various terms by 
using toy models which allow understanding the final results. 

Therefore, the aim of this article is in part to understand in detail the origin of 
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this result and above all to see whether this approximation will remain warranted in 
the post-Planck era, or to determine at which precision it might break down. To do so, 
we compute analytically the lens effect using a perturbative expansion in terms of the 
magnitude of the deflection field and give numerical results for the dominant order. 



2. Weak lensing of the Cosmic Microwave Background anisotropics 

CMB photons, as they emerge from the last scattering surface, are subject to the 
weak gravitational lensing effect. Measurements of the temperature anisotropies, 
and polarization patterns are perturbed by the cumulative effects of the large scale 
structure gravitational wells from z ~ 1000. The net result is that our measurement of 
temperature or polarization of the CMB in a direction n in fact provides an information 
on photons emerging from the last scattering surface in the direction n + £, £ being the 
deflection induced by the gravitational shear 

f(n)=T(n + 0, (1) 
Q(n) = Q(n + 0, (2) 
U(n) = U(n + ®. (3) 

In this paper, we denote by A the lensed appearance of any field A. The deflection £ 
is the cumulative effect of all deflections produced by each gravitational well crossed 
during the propagation of the photon. Since those deflection are small compared to the 
size of the CMB anisotropies, it is sufficient to compute this cumulative effect on the 
unperturbed path of the photon. This approximation is often referred to as the Born 
approximation. By doing so, we ignore higher order corrections to the weak lensing 
effect, such as extra non gaussian correction of the lensing field due to the lens-lens 
coupling and apparition of a curl component in the spin-2 shear field. Both have been 
shown to be negligible in the CMB context [T9] . 

In the Born approximation, the deflection field reduces to [10] 

£(n) = -2 f d X / y / Xcm *7f\ w(xn; r, - x ) , (4) 

JO JMXcmbJjKU) 

where fic(x) is the comoving angular diameter distance and \I> the gravitational 
potential, linked to the density perturbation through the Poisson equation [3]. 

The deflection £ being small, one can evaluate accurately the effect on the CMB 
field by using a perturbative expansion, 

X(n) = X(n) + &3*X(n) + ~^c>W(n) + O^ 3 ) , (5) 

where X is either the temperature anisotropy field T or the polarization components Q 
or [/ pi [20]. 

In the following, we will be interested in the effect of the gravitational shear on 
the properties of the power spectra of the CMB. Since our goal is only to evaluate 
the contribution of the non-gaussian corrections due to weak lensing, and since these 
corrections arise mostly at small scale, it is sufficient to work in the flat space 
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approximation and use a decomposition in Fourier modes rather than in spherical 
harmonics Y^. Let us introduce the lensing potential (f>(n) such that £ = V0(n) and its 
Fourier transform 

0(1) = J dn0(n)e~ a - n . (6) 

([SJ) can then be rewritten as 

X, = Xl 0) + X® + X[ 2) + 0(e) ■ (7) 

When X is a scalar quantity (namely for temperature), the first and second order terms 
now read [20] 

= - J xao 0(1 - ii) [ii • (i - ii)] , (8) 
^ 2) = \ S ^r^w^wa-ii-Wi-ypii-a-ii-ia)]^) 

Finally, it is convenient to use another description of the polarization field than 
the usual Q and U Stokes variables. Indeed, these directly observable quantities have 
awkward geometrical properties. They transform under rotation as the components of 
a spin-2 vector. One then introduces a decomposition in terms of a gradient and curl 
component, named E and B in analogy with the electromagnetic field decomposition, 
which have simpler properties under angular transformations. E is a scalar field and 
B a pseudo scalar one (meaning that the B polarization changes sign under a parity 
transformations). This decomposition has the advantage of offering an important test 
of the primordial cosmology model. The polarization sky pattern essentially maps the 
local quadrupolar temperature anisotropies on the last scattering surface and at the 
linear order, the scalar metric perturbations can only produce an E polarization [21 J . 
The tensorial contribution from gravitational waves is therefore the unique source of 
primordial B polarization, which is concentrated at large scales. Since the tensorial 
contribution yields only a weak contribution to the temperature and E-type anisotropies, 
which is therefore very hard to detect or differentiate from the scalar part, the B modes 
therefore offer a potentially unique opportunity to detect a background of primordial 
gravitational waves, at least on the largest scales. The lensing induced B polarization 
may nevertheless hide the primordial contribution and is in any case dominant at small 
scales [nag. 

The E and B decomposition in the flat sky approximation is easily computed from 
the Stokes variables through [20] 

± X{n) = Q(n)+iU(n) . (10) 

Since these ±X are spin-2 quantities, their Fourier coefficients contain an additional 
factor exp[±2z(<£>/ 1 — (pi)]. We recover the E and B Fourier coefficients with 

±X(l) = E(l)±iB(l) , (11) 
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with which we obtain that the Eqs. ([312]) translate into 

X0(1 -IJih - (1 - li)] , (12) 



= - / WO cos ( 2< ^) + WO M^Ol 

X0(1 - l x ) [l x • (1 - l x )] , (13) 
for the first order and 

^i (2) = 5 // d j0 L WO oob(2 W ) - S(10 sin(2^)] 

x0(l 2 )0(l - li - 1 2 )(1 1 • 1 2 ) [lx ■ (1 - h - 1 2 )] , (14) 



5i (2) = ^ // WO cos(2^0 + E(h) sin(2^)] 

x0(l 2 )0(l - h - 1 2 )(U • 1 2 ) [l x • (1 - lx - l a )] , (15) 

for the second order in the perturbative development. We have used the following 
definition for the angles, 

(px = <P h - (pi ■ (16) 

We can now turn to the computation of the power spectra of the CMB, C* Y , 
defined by 

(X*Yir) — (27r) 2 6(1 — 1') Cf Y , (17) 

where X and Y can be either T, E or B. Since no confusion can arise, from now on 
we shall drop the tilde on the lensed field in order to lighten notations. Putting (flTf) 
and ([HHHD in ([ZD, applying Wick theorem, and truncating the results to second order in 
(f) gives [20] 

cf T = <r (i - yi) +j^J crqf • (i - ii)] 2 d 2 ii , (is) 



where we have defined 



POO 

Jo 



(19) 



With the same computations, the polarization and temperature-polarization cross 
spectra reads 



Cf E = Cf E (l - ^ +j^y 2 J d 2 li [CfW^) + sin 2 (2^ 

xC*t hl [h-(l-h)] 2 , (20) 
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Cf 3 = CP B (l - ^ +j^sj d 2 li PfW(2^) +Cf* sin 2 (2^)] 

xc^tu-a-ioi 2 , (21) 

and 

Cr = Cr(l-^i) +j^y 2 J 008(2^)0^- (l -h)?d\ . (22) 

3. Covariance matrix of the power spectra 

3.1. Gaussian assumption 

Lensing is quite well known to produce non-gaussian features [22j [23J, [21]. Still, in 
previous literature, the weak lensing effect has usually been taken into account in the 
partial way we recalled above. The lensing corrections have been applied to the power 
spectrums only; no effort have been made to reproduce the deviation from gaussianity. 
In other words, the effect of lensing has been reduced to a modification of the damping 
tails of the temperature and E power spectra, and the apparition of the small scale B 
induced by lensing. In that case the power spectrum covariance matrix is diagonal and 
reduces to the square of the Cgs of type TT, TE, EE and BB . 

This approximation might turn out to be a valid one. This can be the case if the 
deviation from Gaussian behavior induced by lensing is small compared to the dominant 
Gaussian contribution, at the level of accuracy of the planned experiments. A. Cooray 
has argued qualitatively that this should be the case for temperature anisotropics [25J. 
A. Lewis has shown numerically with some mock data that this approximation holds 
when constraining the cosmological model at the level of precision of the upcoming 
Planck experiment [16J. The latter approach has the great merit that it does not make 
assumptions on the accuracy of the perturbative expansion used to compute the weak 
lensing effect on the CMB. Moreover, one can easily incorporate non-linear evolution of 
the matter density in the simulation and produce very accurate predictions. A downside 
is that such method has a relatively high computational cost and is physically less 
transparent. 

We propose here another approach to validate the simplification described above. 
We evaluate analytically the impact of the non-gaussian component of the covariance 
matrices of the power spectrum. This provides us with a direct estimate of the extent to 
which the diagonal covariance matrix approximation is a valid one. To do so, we do not 
compute the full non-gaussian correction, and restrict ourselves to the dominant order 
of the lens effect within the perturbative framework introduced in the previous section. 

In the simplest cases (e.g. full sky, noiseless case), one can simply estimate the 
cross-correlation power spectrum of two fields X and Y, by 

d ? Y = w T, ^ + x i r >) = v £ • (23) 
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since the Qs are real. Here X,Y = T, E, B, |1| = £ and Vi is the volume of the sample 
of scales I accessible by the experiment. These estimators of the power spectra are 
unbiased (in the absence of noise), and therefore (cf Y ^ = Cf Y . 

The covariance matrix of the power spectra describes how two estimators of the 
power spectra are correlated. It is given by 

Cov(£, t) XX _yy = (C^CV) - (d? x ) (6F) . (24) 

In the case of a Gaussian random field X, and for a noise free experiment, the Wick 
theorem gives 

Cov(£, nxx^xx = £ [Cf x ] 2 6(£ - £') , (25) 

V£ 

In the case where the field is sampled on the full sky, the volume Vi reduces to 
21 + 1. Neglecting the exact shape of a survey^], this function is often approximated by 
(2£ + 1)/ S ky, the / s ky factor being the ratio of the survey area to the full sky. 

When one corrects for the experimental beam Be, and taking into account 
experimental noise, modeled by non-correlated white noises, one reproduce the usual 
formulae for the power spectrum covariances [271 123 29J 

where X = T, E, B and w^ p denote the power spectra of non-polarized and polarized 
noise respectively. The diagonal term for TE is 

Cov(£)™ = (2£ + 2 1)/gky { {Cj E f + [CJ T + (wrB^} [Cf E + {w P B^] } . (27) 
Finally, there are three non-diagonal terms which are non zero 

c °^ TT - EE =w^ y CIE2 > (28) 

Co,(ef T - TE = (2e + \ )f Cj E [CT + KB,)" 1 ] , (29) 

and 

Cov(t) EE -™ = {2£ + \ )fsk Cf E [Cf E + M)" 1 ] • (30) 

§ Of course, this is only an approximation valid in the limit of very small scales (large £) as compared 
to the size of the survey. When this is not the case, different Ps can be correlated, see the detailed 
work of l26l. 
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3.2. In presence of lensing 

3.2.1. Temperature If we take lensing into account, the computation above is no more 
valid. Indeed, we can no longer assume that the temperature and polarization obey 
to Gaussian distribution. The first term in the r.h.s. of (|2"%|) will translate into a four 
point moments that does not reduce to the usual Gaussian case through Wick expansion 



Of course, we don't have to deal with all four point functions, but only with the 
reduced quantity 

(XiTfUvVf) - (Xtf) {U V V X 1) . (31) 

We use the same perturbation approach than in the previous sections, and reduce our 
computations to second order in lensing for all polarizations but BB — BB. This latter 
term is a particular case, for which the expansion need to be done until order four as 
we will develop later. We do not expect higher order terms to modify significantly our 
results. 

To calculate at order 2 in lensing the covariance for all possible polarization, we 
replace each X, Y, U and V in the formula above by their second order lensed version 
for XY and UV taking the values TT, EE, BB, TE. Then assuming the unlensed 
temperature anisotropies and polarization are Gaussian, one can apply Wick theorem 
to compute the covariance. To simplify the computation, and since we expect the most 
important contribution to arise at small scale, we will keep the flat sky approximation, 
and assume that the volume of sample is the full plane. The complete development is 



given in Appendix A We only summarize the result here. 

For all XY = UV covariances, we obtain two different terms. One being non-null 
only at £ — £' that we refer in the following as the diagonal or Gaussian term. The other 
term, being non-zero for i ^ £' will be called non-diagonal or non-gaussian. In the case 
of the TT — TT covariance, those reads respectively 

V TT-TT = 2 ^ _ ^ + _±_ C TT J d% C^Cf^ih ■ (1 - U)] 2 , (32) 

Ki>~ TT = 2 J ^ C t-v\ {CF [!' ■ (1 - 1')] + Cr [1 ■ (1' - 1)] } 2 • (33) 

We recognize in (1321) something very similar to the classical result. The 
diagonal part is simply the square of the lensed power spectrum, fflHl) . truncated 
to second order in lensing. This is the approximation commonly used, where the 
power spectrum covariance with lensing are simply computed by replacing the power 
spectra by their lensing counterparts in ( J26ll30l) . This approximation works as if the 
temperature anisotropies and polarization remained Gaussian, ignoring the non-diagonal 
contribution, fl33l . that describes the apparition of non-gaussian features. 
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3.2.2. Effects of experimental limitations We have ignored in (I32H33I) the experimental 
limitations of the considered experiment, such as the sensitivity, resolution or sky 
coverage. Since we are interested in the comparison between the diagonal and non- 
diagonal parts of the covariance matrix, these experimental limitations shouldn't affect 
the present work. Let us show how these quantities will enter the results given in the 
present work. The consequence of the introduction of the sensitivity and resolution is 
that the observed fluctuations are now given by 

I, obs = X X B) 12 + n{l) , (34) 

where X = {T, E, B}. We have assumed a Gaussian, uncorrelated noise n(l) with power 
spectrum denoted 

W'Kffll^P-iV 1 (35) 

and the beam transform B[ is related to the FWHM of the beam 9^ through 

B l = exp [-1(1 + 1)02/8 In 2] . (36) 
As a consequence, 

(x obs (x obs r) = 6{l _ y) (p iBi + u -l (/) j ? (37) 

and the estimator of the lensed spectrum reads 

>-l | f -c^obs/ \robs\* „„-l 



a = Bi 1 u -ztx^(x? b r - wr 1 ) . m 

In all the paper we will assume that the temperature noise n^(/) is uncorrelated with 
the signal (and with the polarized noise). Using the definition of the covariance matrix 

Cov(Z,0 = (Wi')-(^i)<&) (39) 
and (|38p . even in presence of the noise, the covariance matrix reads 

Cov(/,/') = (y (x° hs (x° bs )*x^ s (x^ s y) 

- (X^iX^Y) (X£ bs (X£ bs )*)^ Bf 1 B v ~ l . (40) 

We can re-introduce (134]) . and using the fact that the noise is Gaussian and uncorrelated, 
we show that 

Cov(/,/')= [^^(x&tXyX}) +AC l wi 1 B- 1 + 2wi 2 B- 2 , (41) 

J 27T 2lX \ / c 

where 

'xxfXyx;) = (xx;x v x*) - (xx*\ (x v x;) . (42) 



Note that since the lensed anisotropies are not Gaussian, we cannot use the Wick 
theorem to expand the first term of (jHj). Rigorously, it has two contributions one 
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that is diagonal T>i and the non-diagonal term Nip, given by 
we obtain 

Cov(Z, I') = 5(1 - I') [v e + 4Ct (wtBt)- 1 + 2 (tu,S,)- 2 
~ 25(1 -I') \d l + (w l B l )- 1 ] 2 +M,v ■ 



As a conclusion, 



+ Mm 



(43) 



This generalizes (1261) . Note that the previous formulae are easily transposed to 
other polarization terms of the covariance by introducing a cross-correlated spectrum 

estimator C l . We obtain generalizations of ( |26ll30l) . each time introducing the 
quantities T>i and Niy. In the rest of the paper these quantities are calculated and 
compared for every polarization terms of the covariance matrix. 



3.2.3. Numerical evaluation of temperature To evaluate and understand better the 
structure of the covariance matrix, instead of showing comparisons between numerical 
results for ( 13211331) for a few sets of cosmological parameters, we rather show results for 
some simple approximations of the power spectrum, with increasing complexity. 

As a first step, figure [1] shows the comparison when Cg = 5(£ — £ c ), and £ c = 
(200, 400, 800, 1600). This illustrates the effect of lensing on the covariance matrix for 
a single mode. The non-diagonal terms are at least a factor 10~ 4 below the diagonal 
part. We see that the lensing effect spreads the covariance matrix around the l c mode 
in a symmetric way; the amplitude of the effect grows with £ c , the coupling between 
modes being more important at small scales. 



(l c ,r) with C, T = 1 *dirac(l c — 
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2000 



■TT—TT for a 



Figure 1. Non-diagonal (/ ^ I') contribution to the covariance matrix M l l 
Dirac temperature power spectrum Cf T = S(l c — I) and for various values of l c . From 
left to right, l c = 200, 400, 800, and 1600. 



Of course, the temperature power spectra is more complicated than a Dirac 
function. In fact, ignoring the acoustic peaks, it is rather well described by Cj T ~ 
£~ 2 exp ^— ^ V We show in figure [2] the non-diagonal term for this approximation. 
The figure shows 1(1 + 1)1' (I' + l)J\fw in order to see further than the dominant variation 
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Figure 2. Non-diagonal (I / I') contribution to the covariance matrix M l 
temperature power spectrum approximated by C, 



■TT-TT 



for a 



TT 



exp 



values of l c . From left to right, l c = 200, 400, 800, and 1600. 



(-ft) 



and for various 



due to the £~ 2 behavior of the power spectra, as is usually done with the Ce- The figure 
exhibits the same features than the previous one. Namely it spreads over a broad range 
of £, but with a very small amplitude. This spread is more important at large scales 
(small C.) than at small scales, which seems to contradict the idea that weak lensing is 
essentially a small scale effect. Of course, this is due to the exponential suppression of 
the small scales in the power spectrum. We showed above that the effect is essentially 
symmetrical. 

Full results for a concordance model power spectrum are showed in figure IB21 page 
1251 first column. It is also illustrated on figure [3] below. They display the same general 
features we demonstrated on our two simple examples. The structure of acoustic peaks 
(and especially the first one) complicates somewhat the results. However, the fact that 
the non-diagonal contribution is far below the diagonal part remains true. 



3.2.4- Polarization terms We can perform a similar analysis for the i?-type 



polarization. The diagonal and non-diagonal terms read (see Appendix A for a complete 
development) 



V 



EE-EE 



2(0 2 (1-^ 



-a, 



+ 



4 



2tt"V ' (2tt 

BB -2 



EE 



+ c,fw(2^) [lx-a-ix 



(44) 



Af Ej 



EE-EE 



2 / % cosW){cr [l' • (i - 1')] + cr [1 ■ (1' - 1)]} 2 . 



(45) 



Before giving the full numerical result for a concordance model power spectrum, we 
again demonstrate the features of the non-diagonal term on the simplified model 
Cf E = exp ( — ^ ) • figure 0] shows the corresponding result. 
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N TT ~ TT (200, 
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Figure 3. Non-diagonal contribution to the covariance matrix A/^;, TT , normalized 
by the diagonal value £>200~ i an d f° r the temperature power spectrum of the 
concordance model. The cosmological parameters are those measured by WMAP 3 

EH- 



4.0x10 




2000 



EE-EE for 



Figure 4. Non-diagonal (I ^ I') contribution to the covariance matrix Af, , 
a E polarization power spectrum approximated by Cf E — exp ^— J^J and for various 
values of l c . From left to right, l c = 200, 400, 800, and 1600. 



The results for the full concordance model are represented in figure \B2\ second 
column. The main idea is illustrated in figure [5j where one can see the effect of the 
oscillations of the E spectrum and the sub-dominant amplitude of the non-diagonal 
contribution. 

The terms for the B polarization at second order in lensing are obtained by replacing 
E by B and B by E in (14411451) . Surprisingly, there seems to be no non-diagonal 
contributions from the E polarization to the covariance matrix of the B polarization. 
In fact, this is an artifact of the truncation of the perturbative development, which 
explains why J\f l B l P~ BB is so small compared to the Gaussian part (see figure IB21 third 
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N EE ~ EE (200,r)/D(200) 

2.0x10 6 1 




1500 2000 



Figure 5. Non-diagonal contribution to the covariance matrix A/" 200 ,/ , normalized 
by the diagonal value X^jo ~ E ■> and for the power spectrum of the concordance model. 
The cosmological parameters are those measured by WMAP 3 pQ. 



column). Indeed, contributions from the E mode to the B covariance cannot appear 
at second order (in lensing), they appear only at fourth order. Moreover, if we restrict 
ourselves to order 2 for the diagonal parts, there are no contribution proportional to 
(Cf E ) 2 , which may be dominant. 

It is at first sight surprising that, while doing a perturbative expansion, a 4th order 
term may end up being larger than a 2nd order term. The problem comes from the fact 
that this development is in fact done on the Q and U Stokes parameters, where second 
order terms are indeed greater than 4th order ones. However, we are combining these 
Q and U terms to form the E and B fields in a way that enhances the 4th order terms 
relatively to the 2nd order ones in B. The most trivial of those fourth order terms is 
given by the square of the second order lensed B power spectrum which contributes to 
the diagonal part of the covariance 



V 



BB-BB 



(cry 



l 2 2 

2 - -al 



7T 



+ 



(2tt 
2 



n BB 



d^qTui C,r cos^i) + C£* sin^) [U • (1 - U)] 



(27T) 



d\C;« C™ cos^) + C£* sin^(2^) [h ■ (1 - h)Y 46) 



This term, in addition to be the most trivial one, is also the dominant 4th order 
diagonal contribution. Indeed it is a configuration "1+1+1+1", and involves therefore 
a term containing only E modes (contrary to configurations "2+2+0+0", "3+1+0+0", 
or "4+0+0+0"). Note that, as detailed in Appendix A , a configuration is called 
"1+1+1+1" when the four-point correlation function involves four fields at order 1 
in lensing, see e.g. the configurations represented in figure IBll 

We may think that the new diagonal term dominates the diagonal part of the 
covariance matrix, as lower order terms depend on the B power spectra which is much 
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smaller than the E modes. However, on the other hand, this term is of order four in 
lensing and this effect suppresses the term: it finally represents a correction of order 
0.01% to the order 2. 

Let us now turn to the order four non-diagonal contributions. If we assume, 
temporarily, that the B modes are negligible compared to E modes at the same order 
in lensing, we can see that the dominant non-diagonal 4th order terms involve four first 
order E modes from B^\ All possible terms of the form "1+1+1+1" are represented 
graphically in figure IBll 

K B ti B = j^yM + B + c )i ( 47 ) 

where the three terms are given by 

A= . <v 

B = 

and 

C = 



J d\ (O a Cp u| Cp_ u| sin a (2^)8in a (2^) 

x [h ■ (1 - h)f [h ■ (Y - h)] 2 , (48) 
^ f d 2 l!d 2 l 2 Gl E Gl E Cf_ yx ff_ H sin(2^) sin(2^) sin(2y. 2 ) sin(2^) 



x [h • (1 - l x )] [li • (1 - l a )] [1 2 ■ (1 - la)] [1 2 ■ (1 - li)] 5(1 + 1' + l x - l a ) (49) 



d<p' 



J dVl.W (ap li| ) 2 sin 2 (2^ 1 )sin 2 (2^) 

X [lj • (1 - l x )] 2 [1 3 ■ (1 - U)] 2 5(1 + 1' - l! - 1 3 ) . 



(50) 

In these formulae, 5 is the Dirac distribution and we have used the following definitions 
for the angles 



<Pi = <Ph ~<Pi, Vi = <PU ~ W 
Our results agree with those of [301 EH]- 
The TE-TE contributions read, 



(51) 



V 



TE-TE 



+^^V d2i ^r^i ll [ii-a-i 1 )] 2 

+ ^Q Tr /d 2 l 1 [l 1 -(l-l 1 )] 2 [^cos 2 (2^ 1 ) + Qf B sin 2 (2^ 1 )]qr il| 



-pr J axe™ cos^ocpy [u-a-ix) 

for the diagonal part whereas the non-diagonal part is given by 



(27T) 



(52) 



E-TE 



(cr) 2 +c? E cr\ [I'-a-i')] 2 

(cr) 2 + crcr] [i ■ a - n 2 }^ <x*m 
+ [r ■ (i - io] [i ■ (i - 1')] {crcr [i + cos 2 (2^ 



+ 



+ [crcr + crcr] „ 



(53) 
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As already mentioned, until order 2 in lensing, the covariance terms involving 
polarization TB and EB are null. At this order, there are also no cross- correlations 
involving BB polarization. Thus there are no other terms in the diagonal. 

We are left with six off-diagonal (UV 7^ XY) terms. Their analytical expressions 



are given in Appendix B We would like to point out that among these 6 off-diagonal 



terms three of them represent corrections to Gaussian terms ~ Y . But the three oth- 
ers involve polarization of the form BB — XY with XY e {TT, EE, TE}, polarizations 
for which the Gaussian terms were null. This introduces another lensing modification 
to the covariance matrix. 



3.2.5. Numerical evaluation of the corrections Let's now turn to the numerical 
calculation of the non-gaussian terms (ftf l x l Y ~ uv ) compared when possible to the 
Gaussian ones (Vf Y ~ uv ). This is represented in figure lB2llB3l and lB4l These numerical 
calculations mainly extend to all polarizations and to different values of multipoles the 
results represented for TT-TT and EE-EE in figures [3] and EJ On the figures lB2l and lB3l 
we can see that the non-gaussian contributions are completely sub-dominant compared 
to the Gaussian contributions, namely of the order of 10 _2 % or lower. At the end of 



Appendix B in table IBll are given the absolute values of the Gaussian contribution 
to the covariance matrix, in order to calculate the absolute values of the non-gaussian 
corrections represented on figure IB21 and IB31 The term BBBB also receives important 
non-gaussian contributions from the order 4 in lensing. It has been shown in [301 EI] 
that these corrections are negligible at low multipoles but become important at higher 
multipoles (£ > 800). Despite the fact that these terms are of order four in lensing, the 
fact that they involve the E spectrum and not the B spectrum dominates and make these 
terms dominant. This is confirmed by the semi-analytical approach of A. Lewis [16J. 
The lensing also introduce new correlations in the covariance matrix. These terms are 
represented on figure IB4I The amplitude of these new terms is clearly sub-dominant 
compared to any Gaussian term in the matrix. 

Our results fully agree with those of [T71 [18] , although they both assumed that the 
primordial B modes are vanishing. However, in our calculations, we found out that this 
approximation is not recommended in the sense that this arbitrarily sets the second order 
non-diagonal contributions to the BB-BB covariance to zero. For low multipoles, these 
terms represent the leading corrections to the Gaussian assumption. On the other hand, 
one should turn to [T7J [18] do see how the non-gaussian corrections propagate to the 
errors on cosmological parameters, and the consequences on the observation strategies. 

Conclusions 



CMB anisotropies and polarization data are a powerful tool to constrain the cosmological 
model. Indeed, their statistical properties can be computed and compared to the actual 
data. In the minimal case, when weak lensing is neglected, the theory predicts that, at 
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dominant order, the anisotropies and polarization must obey a Gaussian distribution, 
thus allowing the well known data compression that reduces the experimental data to a 
set of power spectra. This is why most article put a strong emphasis on the evaluation of 
the power spectra of the CMB, taking into account or not secondary effects, and reduce 
all experimental results to a set of C/s. This can turn out to be a poor approximation; 
by doing so, one would ignore any deviation from the Gaussian behavior that can arise 
from those secondary effects. 

We have computed analytically and numerically the Gaussian part of the covariance 
matrix as well as the non-gaussian contributions due to lensing. This last contribution 
was usually assumed to be negligible. We prove that this assumption is justified, for 
all polarizations except BB-BB and independently of the sensitivity of the experiment 
considered. The error made is always completely sub-dominant, of the order of 0.01% or 
lower. The covariance matrix can thus be computed using the Gaussian assumption, as 
described in sec. 13.11 : the covariance matrix can be computed by assuming that lensed 
CVs are Gaussian and by using the unlensed formulae (12611301) . but for the replacement 
of the CVs by the lensed spectra C;'s. We can see with [43] that our conclusions remain 
valid independently of the considered experiment, even the most sensitive one. These 
results confirm the recently published [TTl [T8] where the lensing effect on covariance has 
been studied as well. 

The case of BB-BB polarization requires a more extended expansion in lensing, until 
order four, in order to take into account all dominant effects in non-gaussian corrections. 
Indeed, at order four, terms where only E modes contribute to the covariance matrix 
can exist. They have been found [30] (see also [31]) numerically dominant over Gaussian 
terms only when considering multipoles higher than i > 800. 

The weak lensing has another effect on the covariance matrix. It introduces new 
correlations between BB and TT, EE, and TE spectra. Their amplitude is strongly 
sub-dominant compared to any Gaussian contribution to the covariance matrix. 
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Appendix A. Graphical representation of different contributions to the 
four-point correlation function 

The covariance matrix is proportional to the connected part of the four-point correlation 
function of the lensed anisotropy fields, defined b}{jj 

(x l Y l *U ll V^) c = (xftUvVj) - (xtf) (u v V?) . (A.l) 

We have chosen the same Fourier transform conventions as W. Hu [201 122], in which, 

(X h Y*) = (2tt) 2 5(/ 1 - l 2 )C* Y = (2n) 2 5(h - h)C* Y , (A.2) 

for X, Y = <f>, T, E, B and the four-point correlation function may be written in order 
to define a diagonal x>f Y and a non-diagonal part Hf Y ~ vy as 

(xffUyV;} = (2vr) 4 5(/ - l')Vf Y - uv + (2ir) 2 Af l x /- uv . (A.3) 

For simplicity, we will illustrate our method of calculation on the XY — UV = 
TT — TT polarization term in the covariance matrix which will clarify the various 
contributions to the four-point correlation function. The generalization to other 
polarization terms is straightforward. The method we describe here can be extended to 
all orders in lensing and all n— point correlation functions. 



Appendix A.l. General method 

Listing and calculating all contributions can be done graphically [32]: the task is then 
simpler and more efficient. For this, we can use the analogy with Feynman graphs. 
We will focus on the four-point correlation function of the temperature, up to order 
2 in lensing. Let us represent by a cross x a field T(£) and by a circle o the lensing 
potential <fi. The term we calculate contains two fields T at the multipole 1 and two at 
the multipole 1'. We must remember that in general, each field T} is the sum of lensed 
temperature field lensed at each order 

oo 

7} = 7] + ]TTf , (A.4) 
i=i 

where the first T"W are given in (IHl fToT) . This development can be represented graphically 

as 

fj = x + xo + xoo + (A. 5) 

As a consequence, to consider all contributions we can proceed order by order in 
lensing. Once the order is set (say to n) , all possible combinations of the four fields T are 
written such that the sum of the four orders is equal to n. All ways of correlating them 
(represented by a solid line) give all the graphs contributing to the order n. Considering 
that we calculate the covariance matrix given in flA.lj) . among the fields T, two must 



be at multipole /, two at multipole I' and one field at each multipole carry a * for 
|| This definition follows the one used in previous works, for e.g. [251 fTTl I5U] . 
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the complex conjugate. Then we need to use the the following rules (for temperature) 
derived from fl8l TT5|) . 

x=T, 

xo = - J T(h) 0(1 - 10 [h ■ (1 - 10] , (A.6) 

x ° ° = \ J d ~WF T(ll) 0(12)0(1 " 11 " l2) (l1 ' h) [l1 ' (1 " 11 " la)] ' 

where do is given by (|19p . Finally, we use (1A.2I) to calculate the contribution from the 
graph to the four point correlation function. Two classes of terms will appear : the 
diagonal ones proportional to 5(1 — V) that contribute to x>f Y ~ uv and the others that 
contribute to M{\7~ uv [Note the factors (2%) due to the definition flA.30 ]. To calculate 
the contribution to the covariance matrix, we need to introduce the integration over the 
angles of 1 and 1', 

dip f d<p' 



2tt J 2tt ' 

equivalent to the summation over m and m' in the full sky calculation. 



Appendix A. 2. Temperature at order 

At order 0, the graphs will not involve any o; they are given in figure |A~T1 We arbitrarily 
chose to put fields at multipole I in the left part of the graph and fields at multipole I' 
in the right part. We note that the two graphs on the first row do connect the different 
multipoles: they are called connected graphs. The graph on the second row is said not 
connected and will not contribute to the covariance matrix. 



X x 


X X 


^x 
1 p 


X x 
1 r 


X 

J J 

1 p 






X Temperature field 





Figure Al. Configurations of correlations contributing to the order of the covariance 
matrix. 
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The expression for a graph, for instance the first one of figure IA14 is simply the 
product of the two correlations 

CWXTfT,,) = [(2tt) 2 5(/ - OCT]' • (A.7) 

This graph will contribute to V>i. We can check that the second graph of the first row 
gives the same contribution and thus the total contribution to the covariance matrix 
of graphs of order is twice the previous expression. This yields the very first term in 

Appendix A. 3. Temperature at higher order 

The contributions of order one (as well as all other odd orders) are vanishing since we 
are considering fields T and <ft that are uncorrelated and of vanishing mean value 

(T«j> v ) = (Ti) = (fr) = . (A.8) 

Thus, in graphs with only one o, the lensing field cannot correlate to any other fields 
and one gets a term proportional to (0) which is null. For the same reason, all odd-point 
correlation functions of the CMB are null. 

Let us turn to the order 2 in lensing. Graphically, we need now to add two o in 
all possible configurations (each field T can be expanded until order two). The easiest 
possibilities are to consider three unlensed temperature fields and one order 2 field. We 
will call this first class of graphs, the order "0+0+0+2" . Then, similarly to the figure |A~T| 
there are three ways to correlate these configurations since the two lensing fields must 
correlate together. For a chosen lensed field, the contributions are given in figure IA2I 
The nine other contributions to the same order are identical but with permutations of 
the lensing fields. 

As above, the configuration of the second row does not contribute to the covariance 
matrix. To calculate a term, we simply read which fields are correlated and use ( 1A.6I) . 
For example, the first term reads 

~1 1 



(A.9) 



(Tf ) TS)(ITT l ,)= [(2nf 6(1 - l')Cr] 

x(li-l a ) [li ■ (1' - li - l a )] 

= -\{2,f5{i-i'){cryyi ) (a.io) 

We can equivalently place the two lensing fields on all four temperature fields, and we 
can check that the second graph of first row of figure IA2I give the same contribution. 
Thus, the total contribution to the covariance matrix of graphs of order "0+0+0+2" is 
8 times the previous expression of (1A.10I) . Finally, this term contributes to T>i and gives 
the second (negative) term of fl32|) . 



Still restricting to the order 2 in lensing, we are now left with all the possibilities 
involving two temperature fields of order 1 and two unlensed fields. These graphs will 
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Figure A2. First set of configurations of correlations contributing to the order 2 
involving three unlensed temperature fields. More precisely, they contribute to the 
order "0+0+0+2". 

be called the order "0+0+1+1". To obtain all possible terms, we write graphs for all 
possible way of placing the lensing fields o : there are 6 possibilities, represented on 
figure IA3I 

Then for each graph, we write all possible ways of correlating the fields: there are 
still three possible ways. In total, there are 18 graphs of order "0+0+1+1", two of them 
being graphs that do not contribute to the covariance matrix (unconnected graphs). 
Four of them are "not topologically connected" which means that they correlate the 
two multipoles but can be separated into two different two-point correlation functions: 
these graphs then contribute to T>J T ~ TT . They give the third and last term of (1321) . 

We are now left with 12 graphs that are either "connected" or "topologically 
connected". These contributions will give new graphs that cannot be obtained via 
the product of lensed Cis. Four of them are vanishing. The eight other graphs are 
represented in figure IA4I 

The two graphs of the first (resp. second) row give the contributions to MJJ~ TT 
proportional to (Cf T ) 2 (resp. (C Z T T ) 2 ) in ( 1331 ). The last two rows give the contributions 
proportional to Cf T Cf T . 



Appendix A. 4- Rules for the polarization 

Let us finish this section with the few explanation for the polarization case. The method 
explained at the beginning of Appendix A are still valid except that the rules are slightly 
modified. If + (resp. — ) denotes the primordial E (resp. B) field, the rules for 
temperature are replaced by 

El = +E + +E° + ~E° + +E°° + ~E°° + ■■■ , (A.ll) 
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Figure A3. All possible configurations of temperature and lensing fields that 
contribute to the order "0+0+1+1". 



where for example +e° (resp. —e°) denotes the E (resp. B) modes contributing to E. 
Once the graphs are written, their contributions can be calculated using the following 
rules, derived from (11211150 

+e = E , 

;2 



+EO = ~ I E{h) cos(2v?l) 0(1 " h) [h ■ (1 ~ ' 



<N 1 B(h) sin(2^) 0(1 - l x ) [l a ■ (1 - h)} , (A.12) 



(2vr) 2 

+ E oo= l -j E(h) cos(2^) 0(1 2 )0(1 - l a - l 2 )(h ■ 1 2 )[1 1 • (1 - h - h)] , 

E o o = — J °LJ_JL B{h) sin (2^) 0(1 2 )0(1 - 1, - h )( h . l 2 )[ h . (1 - l x - 1 2 )] . 

Using these rules, all the procedure is strictly identical to the temperature case. For 
example, the graphs contributing to the polarized terms of the covariance are given by 
the same diagrams given in figures IA11 IA2I and IA4I The only difference is that each 
field contain now 2 terms from the E modes and the B modes. However a lot of terms 
are vanishing due to the fact that 

(T h BQ = (E h BQ = . (A.13) 
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Figure A4. Non vanishing, "topologically connected" , configurations of correlations 
contributing to the order "0+0+f+l". These graphs are all the contributions to 

K r ~ TT - 



We can write the same expressions for the expansion of Bi, by replacing E by B and B 
by E 

Bi = -b + -bo + + B o + -goo + + B oo + ••• , (A.14) 
At every order, clearly the main contribution will come from the the E modes denoted 



+b°, +b ° ° • • •• This diagrammatic representation is employed in section Appendix B.2 



for the calculation of contributions to the covariance BB-BB at order 4 in lensing. 



Appendix B. Full expressions for the covariance matrix 

Appendix B.l. Cross-correlated terms of the covariance matrix at order two 

We give here the expressions of the off-diagonal terms of the covariance matrix. Non- 
gaussian contributions A^[, y-X are numerically evaluated and represented on figure IB~3l 
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They are compared to the Gaussian part V^ V ~ XY , which is usually taken into account. 

r j?B =2(o 2 (i-^) 

+ Wy CfE I d2l i^ E cos(2 (/?1 )C | f ll |[li • (l-li)] 2 , (B.l) 

A/"™ = 2 f%{ ^ [1 ' ' (1 " r)]2 + ^ [1 ' (1 ' " 1)f C ° S2(2V?,) 

+2C™C^cos(2 ¥ /) [V ■ (1 - 1')] [1 • (V - 1)]}CJ* ,, . (B.2) 

V TT-TE = 2 qTT qTE ^ _ + I d2li C TT C H hffli . (1 _ ^2 

+ W? GJT I dX C ™ C0 < 2 ^KU l i • (1 - li)] 2 , (B.3) 

Kr TE = 2 / [i' • (i - 1')] 2 + crcr 008(2^) [i • (i' - 1)] 2 

\crcr + c^cr \v . (i - m [i • <r - i)i . (b. 4 ) 



+ 



The last polarization terms are given by 



^EE-TE r) /~iTE /~iEE ( i ^ Jl\ , 2 f , 2 i \ r<EE ^ c ,2 ( ,„ \ 

V e =2C l C t ^1 - — a Q J + J^nf C i J d L i [ C h cos ( 2 V?i) 



+ Qf B sin 2 (2^ 1 )lc | f ll| [l 1 -(l-l 1 )] 



(2 



^Cf* | d 2 U Cl E cos^C^ [lx • (1- U)] 2 , (B.5) 



M EE-TE = 2 J ^{ C EE C TE CQS 2 M y . (1 _ ^2 

+cr<z E <**w) [i - (i' - 1)] 2 

+ [Cf^CfW^') +C^CP cos(2^)] 

x[]!.{\-X)][\.{)!-\)]}c^. (B.6) 

There are finally terms that are vanishing if we restrict ourselves to Gaussian terms 
only, involving cross terms with the B polarized spectrum. They are given by 

j^BB-TT = 2 J sin 2 MC TE c EE y . (1 _ ^2 (R?) 

^J 8 "** = 2 J ^C** V i{c? B [!' ' (1 " 1')] + Cf E sin(2^) [1 ■ (1' - 1)]}' . (B.8) 
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Kf' TE = 2 / %<t_ n sinW^Cf 5 [1' • (1 - l')] 2 
+C^Cf B [l-(l'-l)][l'-(l-lO]} . 

Appendix B.2. Dominant contributions to BB-BB at order four 



(B.9) 



If we temporarily assume that the B modes are negligible compared to the E modes 
at the same order in lensing, we can see that the dominant non-diagonal 4th order 
terms involve four first order E modes from B^\ As explained in Appendix A , they are 
denoted +b°- All possible terms of the form "1+1+1+1" are represented graphically 
in figure IB11 For each column, the two graphs have the same contribution, respectively 
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Figure Bl. Non-diagonal "topologically connected" configurations of correlations 
contributing to the covariance matrix, for the polarization BB-BB at order 
"1+1+1+1". These graphs are the dominant terms contributing to A/j^; (4) a ^ orc ^ er 
four in lensing. 



given by @BH5DI). 



Appendix B.3. Numerical calculation of the non-gaussian contributions 

We now turn to the numerical calculation of all these terms. They are represented on the 
following figures [B2UB3t and lB41 When possible (figures IB~2| and lB3j) these non-gaussian 
corrections are renormalized by the Gaussian value T>f Y ~ uv so that the amplitude of 
the graph reflects the amplitude of the correction. These renormalization factors are 
given in Table IB11 This table can be used to evaluate the absolute amplitude of non- 
gaussian corrections of figures IB21 and IB3I and can then be compared to the corrections 
of figure | 
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Multipole 


TT-TT 


EE — EE 


BB — BB 


TE-TE 


TT -TE 


TT -EE 


EE-TE 


£ = 200 


2.6 x 10 7 


11 


7.2 


8.4 x 10 3 


5.2 x 10 4 


10 2 


34 


£ = 400 


2.4 x 10 6 


9.2 x 10 2 


1.4 x 10 2 


2.3 x 10 4 


2.0 x 10 4 


1.6 x 10 2 


3.9 x 10 2 


£ = 800 


5.5 x 10 6 


6.7 x 10 3 


4.9 x 10 3 


9.8 x 10 4 


1.3 x 10 5 


2.8 x 10 3 


4.4 x 10 3 


£ = 1600 


4.3 x 10 5 


1.2 x 10 6 


1.2 x 10 6 


3.7 x 10 5 


1.1 x 10 4 


2.8 x 10 2 


1.9 x 10 4 



Table Bl. Absolute values of the covariance matrix in the Gaussian limit T> XY uv 
with the prc-factor [£(£ + l)/2n} 2 . The unit is (fJ,K) A . 
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Figure B2. On the first row, we recall the form of the TT, EE, BB, and TE power 
spectrum in the absence of lensing. The cosmological parameters are given by WMAP 
3 pQ and we assume a tensor/scalar ratio r = 1. Below is the non-diagonal (I ^ V) 
contribution to the covariance matrix Af[ J l Y~ XY for several polarization and various 
values of I. From left to right, UV - XY = TT - TT, EE - EE, BB - BB, and 
TE - TE, and from top to bottom, / = 200, 400, 800, and 1600. The amplitude is 
normalized to by the diagonal value T>V V ~ XY to have the relative contribution. 
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Figure B3. Non-diagonal (I ^ I') contribution to the covariance matrix M^~ XY for 
several polarization and various values of /. From left to right, UV — XY = TT — EE, 
TT - TE, and EE - TE, and from top to bottom, / = 200, 400, 800, and 1600. 
The amplitude is normalized to the diagonal value T>f v ~ XY to have the relative 
contribution. Above the middle column we represent the cross correlated spectrum 
Cj E \ since it is the only spectrum that generates the term NjJ~ EE . The cosmological 
parameters are given by WMAP 3 [1] and we assume a tensor/scalar ratio r = 1. 
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Figure B4. Non-diagonal (I ^ /') contribution to the covariance matrix M^~ XY for 
several polarization and various values of I. From left to right, UV — XY = BB — TT, 
BB - EE, and BB - TE, and from top to bottom, I = 200, 400, 800, and 1600. The 
amplitude is not renormalized by any Gaussian value since these terms are vanishing 
in the Gaussian limit. On the first line are represented the C;s that play the most 
important role in the non-diagonal contribution below. The cosmological parameters 
are given by WMAP 3 pQ and we assume a tensor/scalar ratio r = 1. 
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